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Abstract 

The results of the computer investigation of the sign changes of the 
difference between the number of twin primes ^(x) and the Hardy- 
Littlewood conjecture C2Li 2 (x) are reported. It turns out that 7r 2 (x) — 
C2Li 2 (x) changes the sign at unexpectedly low values of x and for x < 2 42 
there are over 90000 sign changes of this difference. It is conjectured 
that the number of sign changes of ^(x) — C2Li 2 (x) for x G (1,T) is 
given by \ff j log(T). 
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Let 7r(x) be the number of primes less than x and let Li(x) denote the logarithmic 
integral: 

du 



Li(x) = / -— -. (1) 

Ji log(u) 

The Prime Number Theorem tells us that Li(a;)/7r(a;) tends to 1 for x — > oo and 
the available data (see [2]) shows that always Li(x) > tt(x). This last experimental 
observation was the reason of the common belief in the past, that the inequality 
Li(x) > ir(x) is generally valid. However, in 1914 J.E. Littlewood has shown pQ 
(see also [10]) that the difference between the number of primes smaller than x 
and the logarithmic integral up to x infinitely often changes the sign. In 1933 S. 
Skewes [3] assuming the truth of the Riemann hypothesis has estimated that for 

10 34 

sure d(x) = tt(x) — Li(x) changes sign for some xq < 10 10 . The smallest value xq 
such that for the first time vr(xo) > Li(xo) holds is called Skewes number. In 1955 
Skewes jl] has found, without assuming the Riemann hypotheses, that d(x) changes 
sign at some 

xq < exp exp exp exp(7.705) < 10 10 ° 

This enormous bound for xq was reduced by Cohen and Mayhew [5] to xo < io 1q529 ' 7 
without using the Riemann hypothesis. In 1966 Lehman [6J has shown that between 
1.53 x 10 1165 and 1.65 x 10 1165 there are more than 10 500 successive integers x for 
which 7r(x) > Li(x). Following the method of Lehman in 1987 H.J.J, te Riele [7] 
has shown that between 6.62 x 10 370 and 6.69 x 10 370 there are more than 10 180 
successive integers x for which d(x) > 0. The lowest present day known value of the 
Skewes number is around 10 316 , see [8] and [9]. 

The number of sign changes of the difference d(x) = tt(x) — Li(x) for x in a given 
interval (1,T), which is commonly denoted by v(T), see [TO], was treated for the 
first time by A.E. Ingham in 1935 [11] chapter V, [12] and next by S. Knapowski 
[13]. Regarding the number of sign changes of d(x) in the interval (1, T), Knapowski 
[T3] proved 

v{T) > e~ 35 log log log log T (2) 

provided T > exp exp exp exp (35). Further results about v(T) were obtained by J. 
Pintz P3] and J. Kaczorowski [15J 

In this paper I am going to look for the analog of the Skewes number for the 
twin primes. 

Let iT2(x) denote the number of twin primes smaller than x. Then the unproved 
(see however [17]) conjecture B of Hardy and Littlewood [16] on the number of prime 
pairs p,p + d applied to the case d = 2 gives, that 

f x u 

tc 2 {x) ~ C 2 U 2 {x) = C 2 - — 2-—:du, (3) 

J a log [U) 

where C 2 is called "twin constant" and is defined by the following infinite product: 

C aS2 n(l-^I)*) = 1-32032... (4) 
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Usually the lower limit of integration a in ([3]) is chosen 2, but the author believes 
that the proper choice for the lower limit of integration should be 5, not 2, because 
(3,5) is the first twin pair. (Analogously in flTJ the lower limit of integration is 2, to 
ensure that Li(2) = 0.) 

For the first time the conjecture ([3]) was checked computationaly by R. P. Brent 
[18J. This author noticed the sign changes of the difference 7r 2 (x) — C 2 Li 2 (:r) but did 
not elaborate about this further. I have looked on the difference <i 2 (x) = 7r 2 (a;) — 
C 2 Li 2 (x) using the computer for T up to 2 42 « 4.4 x 10 12 . Like for usual primes 
initially C 2 Li 2 (x) > tt 2 (x), but surprisingly, it turns out that there is a lot sign 
changes of d 2 (x) = 7r 2 (x) — C 2 Li 2 (x) for x in the range (1, 2 42 ). For the case of the 
lower integration limit a = 2 in ([3]) the first sign change of d 2 (x) appears at the 
twin pair (1369391, 1369393). However, for the choice a = 5 the first sign change 
of d 2 (x) appears already at the pair (41, 43)! Next "Skewes" pairs for a = 5 are 
(6959, 6961) and (7127, 7129). After that there is a gap in the sign changes up to 
(1353257, 1353259), what is comparable to the case of a = 2. Such a behavior is 
reasonably, because the difference between a = 2 and a = 5 is relevant only at low 
x and the contribution stemming from the interval (2,5) is becoming negligible for 
large x. 

Let v 2 {T) denote, by analogy with usual primes, the number of sign changes of 
d 2 (x) in the interval (1, T). The Table I contains the recorded number of sign changes 
of 7r 2 (x) — C 2 Li 2 (x) up to T = 2 22 , 2 23 , . . . , 2 42 for both choices of the parameter a. 
The values of T searched by the direct checking are of small magnitude from the point 
of view of mathematics, but large for modern computers. The observed numbers 
v 2 (T) behave very erratically, see Fig.l, in particular there are large gaps without any 
change of sign of the d 2 (x). If one assumes the power-like dependence of v 2 {T) then 
the fit by the least square method gives the function aT@, where a = 0.0766741 and 
(3 = 0.479031 . . .. Instead of such accidentally looking parameters for the pure power- 
like behavior (especially a has a very small value) after a few trials I have picked the 
function VT / log(T) as an approximation to v%(T) as a much nicer function without 
any random parameters involved and simultaneously taking values very close to the 
fit aT 13 , see Figure 1. Thus we state the conjecture: 

z^ 2 (T) = VT/ log(T) + error term (5) 

Let us stress in favor of (jSJ) that there are 7 crosses of the curve yT / log(T) with 
the staircase-like plot of v 2 {T) obtained directly from the computer data. The 
last column in the Table 1 contains the values of the function yT/\og(T). If the 
conjecture (JSJ) is true, then there is infinity of twins. Also if (jSJ) is valid, it means 
that the estimation ([3]) is in some sense more accurate than (JTJ), because there are 
more points, where ([3]) exactly reproduces the actual number of twins — for (JTJ) 
there is much less such values of x that Li(x) is equal to ir(x), see (j2J). However, 
presumably the (unknown) error term in ([3]) is larger than error term for ir(x). 

The difference of many hundreds of orders between values of x such that tt(x) — 
Li(x) and 7r 2 (x) — C 2 Li(a;) changes the sign for the first time seems to be very 
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astonishing. Let me give an example from physics: the energy of the ground states 
of the hydrogen and helium are respectively -13.6 eV and -79 eV and do not differ 
by hundreds of orders! 

I have tested the numerical results using several different computers, programs 
and compilers. In particular, to calculate the integral Li2(x) I have used the 8-point 
self-adaptive Newton-Cotes method and the 10-point Gauss method. This integral 
was calculated numerically in successive intervals between consecutive twins and 
added to the previous value. It seems to be natural that different methods gave 
exactly the same values of Li 2 (x) since the integrand in is a very well behaved 
function. At least up to T = 2 32 (this limitation stems from the fact that some 
compilers did not allow larger integers than 2 32 ) all results obtained by different 
runs were exactly the same. 
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TABLE 1 

The number of sign changes of d,2(x). The case a = 2 is in second column, and third 
column contains data for a = 5, while the last column presents values obtained from ([5|). 



T 


^2 CO for a = 2 


v-i (T) for a = 5 


VT 


2 22 


29 


32 


134 


2 23 


29 


32 


182 


2 24 


29 


32 


246 


2 25 


29 


32 


334 


2 26 


238 


269 


455 


2 2V 


854 


942 


619 


2 28 


1226 


1401 


844 


2 29 


1226 


1401 


1153 


2 30 


1226 


1401 


1576 


2 31 


1226 


1401 


2157 


2 32 


2854 


3045 


2955 


2 33 


7383 


7358 


4052 


2 34 


9115 


8974 


5562 


2 35 


12682 


12431 


7641 


2 36 


23634 


23103 


10505 


2 37 


31641 


30770 


14455 


2 38 


31641 


30770 


19905 


2 39 


31641 


30770 


27428 


2 40 


38899 


37904 


37819 


2 41 


55106 


54179 


52180 


2 42 


90355 


89768 


72037 
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Fig.l The plot showing the comparison of the actual values of 2^2 CO for a = 2 and 
a = 5 found by a computer search with the conjecture ([SD- 



